Distillation and purification of symmetric entangled Gaussian states 
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We propose an entanglement distillation and purification scheme for symmetric two-mode entan- 
gled Gaussian states that allows to asymptotically extract a pure entangled Gaussian state from any 
input entangled symmetric Gaussian state. The proposed scheme is a modified and extended version 
of the entanglement distillation protocol originally developed by [Browne et al., Phys. Rev. A 67, 
062320 (2003)]. A key feature of the present protocol is that it utilizes a two-copy de-Gaussification 
procedure that involves a Mach-Zehnder interferometer with single-mode non-Gaussian filters in- 
serted in its two arms. The required non-Gaussian filtering operations can be implemented by 
coherently combining two sequences of single-photon addition and subtraction operations. 

PACS numbers: 03.67.-a, 42.50.Dv 



I. INTRODUCTION 

Quantum communication protocols are very sensitive 
to noise and decoherence that accompanies distribution 
of quantum states over any realistic quantum channel. 
The main source of errors in optical quantum commu- 
nication systems is losses, that e.g. limit the range of 
current point-to-point quantum cryptography to about 
one hundred of kilometers Jj . Significant attention has 
been therefore paid in recent years to finding schemes 
that would allow to suppress the noise and decoherence 
in quantum communication. Of particular interest is a 
faithful distribution of pure entangled quantum states 
between two distant parties Ahce and Bob who can sub- 
sequently use them as resource for some quantum infor- 
mation processing tasks such as quantum teleportation 
or entanglement-based quantum key distribution. Alice 
and Bob can achieve this goal by employing a protocol 
known as entanglement distillation 0,11]. They first dis- 
tribute several copies of an entangled state over the noisy 
channel. From the shared copies of a weakly entangled 
and mixed state they then extract a highly entangled 
pure state by means of local quantum operations and 
classical communication (LOCC). Elementary entangle- 
ment distillation and concentration schemes have been 
successfully demonstrated experimentally for entangled 
two-photon states as well as for two entangled op- 
tical modes 04§|- 

Particularly interesting but also complicated is the dis- 
tillation of so-called continuous- variable (CV) entangle- 
ment where we deal with quantum states belonging to an 
infinite dimensional Hilbert space of field modes. The ex- 
perimentally easily accessible CV Gaussian states, which 
possess Gaussian Wigner function, cannot be distilled by 
Gaussian operations only p^| - [l^ . This means that some 
non-Gaussian operation, such as photon-counting fis'l or 
Kerr nonlinearity ^T^ . [l5| is required. An iterative CV en- 
tanglement distillation protocol for Gaussian states was 
proposed by Browne and coworkers [iTj . Their ap- 
proach combines two steps, see Fig. 1. First, Gaussian 
states are de-Gaussified by an appropriate quantum filter. 
Then, the states are re-Gaussified by an iterative proce- 



dure whose each steps involves interference of two copies 
of a two-mode state on balanced beam splitters, followed 
by projection of one output on a vacuum, c.f. Fig. 1(b). 
It can be rigorously proved that this protocol converges 
to a Gaussian state. This state, however, will generally 
be mixed, only under very special conditions the protocol 
will converge to a pure state. The distillation procedure 
thus increases entanglement of the state, but at the same 
time it reduces purity of the state, as we shall show be- 
low. A possible solution of this problem would be to send 
through the lossy channel a very weakly entangled pure 
state whose entanglement could then be increased by dis- 
tillation [3] or probabilistic noiseless amplification [lol - 
[2^ while preserving high purity of the state. However, 
this approach is not applicable in situations where two 
distant parties share some fixed noisy entangled states 
and want to extract pure entangled state without having 
control over the source of the states. 
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FIG. 1: Entanglement distillation scheme for two-mode Gaus- 
sian states p^. [T7|. (a) De-Gaussification of initial two-mode 
Gaussian state pin by local non-Gaussian operations (quan- 
tum filters) Z. (b) Iterative Gaussification procedure whose 
single step involves local interference of two copies of a state 
on balanced beam splitters BSi and BS2 followed by projec- 
tions of one output mode on each side onto vacuum using 
photodetectors Da and D b , respectively. In the protocol, Al- 
ice and Bob use classical communication to announce success 
or failure of local filtering operations and projections. 
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In this paper we propose a modified and extended 
version of entanglement distillation protocol for sym- 
metric two-mode Gaussian states, which asymptotically 
converges to a pure Gaussian state. This is a crucial 
advantage with respect to previousl y prop osed proto- 
cols for CV entanglement distillation |16l - [l8l | and exper- 
imentally demonstrated CV entanglement concentration 
schemes ^ ^] that cannot fully suppress the effect of 
losses and cannot be used to distill pure entangled Gaus- 
sian state from noisy weakly entangled Gaussian states. 
The GV entanglement distillation scheme considered in 
the present paper has a nested structure. The whole 
distillation procedure including non-Gaussian filtration 
and Gaussification is performed several times, so that 
already distilled Gaussian states are taken as an input 
of the next stage of the protocol. Moreover, we signifi- 
cantly alter the non-Gaussian part of the protocol. The 
de-Gaussification is not performed on a single copy, as in 
previous works [9, 13, 17, 18], but two copies of the state 
are used simultaneously and one of them is consumed by 
the protocol. This makes our approach slightly more re- 
source demanding but also enables us to distill Gaussian 
states with any chosen amount of entanglement and an 
arbitrary high purity. 



II. DISTILLATION OF GAUSSIAN 
ENTANGLEMENT 



Gaussian states can be succinctly characterized by first 
and second moments of quadrature operators. Let Xj 
and pj denote amplitude and phase quadrature oper- 
ators of mode j, respectively. The quadrature opera- 
tors satisfy canonical commutation relations [xj^pk] — 
iSjk- We collect all operators into a vector r — 
{xi,pi,X2,P2, ■ ■ ■ , xn,Pn), where N is the total number 
of modes. The second moments can be collected into 
a a covariance matrix 7 whose elements are defined as 
7jfc = ({Afj, Affe}), where Afj 

notes the anticommutator. The first moments (fj) can be 



set to zero by local coherent displacements which do not 
modify entanglement of the state. In the rest of the pa- 
per we will therefore consider only Gaussian states with 
{r)=0. ^ 

Let us briefly review the iterative entanglement distil- 
lation protocol for Gaussian states proposed by Browne 
et al. |l6l . \vf\ and depicted in Fig. 1. In what follows 
we will focus on a symmetric protocol. More specifically, 
we will assume that a covariance matrix of a two-mode 
Gaussian state p shared initially by Alice and Bob can 
be expressed as 
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The generalized Heisenberg uncertainty relation implies 
that C > 1 and ~ >1 must hold. Furthermore, we 



require that the state is entangled, otherwise any distil- 
lation would be pointless. This yields another constraint 
[l8| , C — \S\ < 1. Pure two-mode squeezed vacuum state 
is recovered for C — cosh(2r) and S = sinh(2r), where r 
is a two- mode squeezing constant. 

As shown by Lund and Ralph [l^ , the entangled state 
([T]) can be conveniently represented as a pure two-mode 
squeezed vacuum state with some squeezing parameter r 
transmitted over a lossy channel with intensity transmit- 
tance T, 



7 = r7TMSV + (1 - T)7vac, 



(2) 



where 7tmsv is a covariance matrix of two-mode 
squeezed vacuum and the covariance matrix of vacuum 
state 7vac is an identity matrix. We immediately obtain 
expressions for C and S in terms of r and T, 



C cosh(2r) + 1 - T, 5 = T sinh(2r) , 



(3) 



whose inversion gives [T^ 



tanh(2r) 



25(C- 1) 
S'2 + (C-1)2 



T = 



2{C-1) 



(4) 

The symmetric entanglement distillation protocol de- 
picted in Fig. 1 preserves the structure of two-mode state 
while altering the parameters C and 5", or, equiva- 
lently the two-mode squeezing r and the channel trans- 
mittance T. As shown in Fig. 1(a), a first step of the pro- 
tocol is probabilistic de-Gaussification of the two-mode 
state by local single-mode filtering operations Z ^ n + w, 
where n is a photon number operator and w is some con- 
stant. Such operation can be implemented by mixing an 
ancilla single photon state |1) with a signal beam on an 
unbalanced beam splitter with carefully chosen transmit- 
tance, followed by projection of ancilla output port onto 
a single-photon state [l3, [1^ • Note that we can write 

Z = (1 — w)a^a -|- wad} , (5) 

and i I de- where a and at denote annihilation and creation opera- 

L'J i- T^T^;^ ^ ^i-^ i-l^^^ „1J- J-* ...U, „ 1 J-U, 



tors. This suggests that, alternatively, one can coherently 
combine two sequences of the elementary operations of 
single-photon addition and single-photon subtraction 
[25l - [27j to implement the desired operation Z [2l,[2l|. A 
feasibility of this latter approach has been demonstrated 
in a recent experiment where a coherent combination of 
operations aa^ and a^d has been implemented in order to 
directly test the commutation relations for bosonic cre- 
ation and annihilation operators [28| . 

After local filtering, the resulting states are then re- 
Gaussified by an iterative protocol. An elementary step 
of this Gaussification procedure consists in local interfer- 
ence of corresponding modes of two copies of the state on 
balanced beam splitters BSi and BS2 followed by projec- 
tion of one output port of each beam splitter onto vacuum 
state, see Fig. 1(b). Mathematically, the Gaussification 
protocol is described by an iterative map [l^, [l3| , 



P 
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Tr [fG(p(") «)/5("))] 



(6) 



3 



where Eg denotes a probabilistic Gaussian operation on 
two copies of a two-mode state p'^^\ Using the labehng 
of modes as in Fig. 1(b) we can write 



£g{p 



ABCD) 



K'^ PabcdK, 



(7) 



where K — J/bs ® ^bs bdI*^' 0)c-D, |0) denotes vac- 
uum state, and Ubs stands for a two-mode unitary op- 
eration corresponding to interference of two modes on a 
balanced beam splitter. If the Gaussification procedure 
converges, then the asymptotic state /5°° is Gaussian, and 
its covariance matrix can be determined from a state /5'^' 
after the first iteration of Eg- We use short hand nota- 
tion Pjk,mn = (j, fc|/5|rn, n) for density matrix elements in 
Fock basis. We define auxiliary matrix 



Poo, 00 



(8) 



Covariance matrix 7out of asymptotic Gaussian state p°° 
is then a function of the matrix elements ajk,mn with 
j + k + m + n = 2. Explicit expressions and details can 
be found in Refs. [TtI. Il8l|. For the class of symmetric 
states ll]) considered here some of the matrix elements 
are equal to zero, 

O'20.00 = O'02,00 = f00.20 = O'00.02 = ClO.Ol = CTOl.lO = 0, 

(9) 

and furthermore iJio.io = croi.oi- As we shall see this 
property is satisfied also by our generalized distillation 
protocol that will be described in Seciton IV. 

A crucial observation made by Lund and Ralph is 
that the distillation protocol shown in Fig. 1 does not 
change the value of a parameter 



e = (1 - T) tanh(r) 



1 



25 



(10) 



It holds that e is the same for an initial noisy Gaussian 
state before distillation pin and the final distilled Gaus- 
sian state p°°. The distillation simultaneously increases 
both effective two-mode squeezing parameter r and ef- 
fective transmittance T while keeping e constant. This 
clearly increases the entanglement of the state. However, 
a simple calculation shows that a purity of the state de- 
fined as V = Tr(p2) is reduced by the distillation. For 
Gaussian states it holds that V = l/y^dct(7). With the 
help of formulas (P), ^ and ^ we obtain 



■P = [l - - cosh(2r) + 2e sinh(2r)] 



(11) 



Since e < 1 by definition, 7^ is a decreasing function 
of r for a fixed e. This shows that the entanglement 
distillation protocol decreases state purity. For symmet- 
ric Gaussian states Q the von Neumann entropy of the 
state depends only on V and is a monotonically decreas- 
ing function of P. So the distillation also increases von 
Neumann entropy of the two-mode Gaussian state. 



III. SINGLE-COPY PRE-PROCESSING 

In order to increase the purity of distilled state while 
preserving (or even increasing) its entanglement, we need 
to find a way how to reduce the value of parameter 
e. We seek to achieve this by a suitable modification 
of the de-Gaussification step of the distillation proto- 
col. We will first consider two simple single-copy strate- 
gies. One option is to replace filter Z = n + w with 
a single-photon subtraction Z = a [l3j . Another op- 
tion is to perform some local Gaussian operations prior 
to de-Gaussification. However, as we shall argue below, 
these approaches do not allow to decrease value of e. In 
Section IV. we therefore develop an alternative collective 
two-copy de-Gaussification scheme that does the job and 
reduces value of e of the final distilled Gaussian state. 



A. Local single-photon subtraction 

If we use single-photon subtraction as the de- 
Gaussification filter Z, then the initial two-mode density 
matrix /5i„ will be transformed as, 



P 



(0) 



ab a^&^ 
Tr[ata6t6pij 



(12) 



where a and b denote annihilation operator of Alice's and 
Bob's mode of shared two-mode state, respectively. For 
symmetric mixed Gaussian state with covariance matrix 
([T]) the matrix elements ajk,mn required for determina- 
tion of final distilled Gaussian state can be analytically 
calculated. The nonzero elements read 



cii.oo — Coo, 11 



ClO.lO — Coi 01 



2TX [1 + 2A2(1 -r)2] 

1-A4(l-T)4 ' 
2T(l-r)A2 [2 + \^{l-Tf] 
1 - A4(l-T)4 



(13) 

where A = tanhr. With these expressions at hand we 
can analytically calculate covariance matrix 7out of the 
distilled Gaussian state p^°°\ The whole protocol pre- 
serves the symmetry of the state and 7out has the form 
([1]), only the value of parameters C and S is changed. 
We are primarily interested in eout- After some algebra 
we find that this parameter can be expressed as a ratio 
of two elements of matrix a. 



Cout 



o'lo.io 
cii.oo 



(14) 



On inserting the explicit formulas ([T^ into Eq. we 
obtain 



tout 



2 + \^{\-Tf 

1 + 2A2(1-T)2' 



(15) 



where €;„ = (1 - r)A. Since A2(1 - T)2 < 1 we have 
Eout > Ein- The replacement of Fock-diagonal filter n-\-w 
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with photon subtraction operation a makes things even 
worse in the sense that the value of e is increased instead 
of decreased. 



B. Local Gaussian operations 

ExperimentaUy simplest filtering would be based on 
application of local Gaussian operations and measure- 
ments. One may conjecture that local Gaussian opera- 
tions might help to reduce value of e although they can- 
not increase entanglement of the state. Let us therefore 
investigate this possibility. Suppose that we perform a 
local Gaussian operation on each mode of the state pin 
prior to the de-Gaussification by filter Z. A general two- 
mode Gaussian operation can be characterized by a four- 
mode covariance matrix T that can be conveniently split 
into output and input parts labeled by indices 1 and 2, 
respectively, 



r = 



Ti ri2 

1^12 



(16) 



This operation transforms covariance matrix of a Gaus- 
sian state according to the formula [Tl| 



y = Ti - ri2 [r2 + S7^bE^] 'r 



12; 



(17) 



where S = diag(l, — 1, 1, — 1). For local operations the 
matrix F reduces to a direct sum of two 4x4 two-mode 
covariance matrices, 



r = F^^, 



BB' 



(18) 



where the unprimed (primed) labels refer to input (out- 
put) modes. We require that the transformation (fT7|) 
preserves the form and symmetry of the covariance ma- 
trix ([l} . This is guaranteed when Taa' = Tb B' and both 
these matrices have the same structure as 7 in Eq. ([!]). 
Moreover it suffices to consider only covariance matrices 
corresponding to pure Gaussian states because any op- 
eration represented by a covariance matrix F of a mixed 
state can be obtained as a Gaussian mixture of oper- 
ations represented by pure-state covariance matrices F. 
This implies that Taa' = ^bb' are covariance matrices 
of a pure two-mode squeezed vacuum state with some 
squeezing constant s. Physically, the resulting opera- 
tion (|17p can be implemented by mixing each mode with 
a vacuum state on an unbalanced beam splitter with a 
properly chosen amplitude transmittance r and project- 
ing one output mode on a vacuum state. The parameters 
C and S' of covariance matrix of a symmetric Gaussian 
state after application of the local Gaussian operations 
can be expressed as follows. 



C 
S' 



C[cosh^(2s) + 1] + (C2 -5*2 + 1) cosh(2s) 



(C-l-cosh(2s))2 
5sinh2(2s) 
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FIG. 2: Two-copy de-Gaussification scheme. BSj denotes 
balanced beam splitter, Z indicates a non-Gaussian operation 
- a Fock state filter Z = h — 1, and each detection block Dj 
makes projection onto a superposition of vacuum and single- 
photon state. 



On inserting these expressions into Eq. ([TU]) we find af- 
ter lengthy but straightforward calculation that e' = e, 
hence this Gaussian operation does not change value of 
e. Another way of proving this statement is based on 
the observation that the considered local Gaussian op- 
eration is a filter diagonal in Fock basis, Zq = t". For 
a symmetric Gaussian state ([IJ the parameter e can be 
expressed as 



Pio.io 
Pii,oo ' 



(19) 



c.f. Eq. (|14p . Density matrix elements of un- normalized 
state after local Gaussian filters Zq are given by p'-^Q = 
^2^10, 10, p'li 00 = t"^ Pii,OQ, and it immediately follows 
that the ratio Pio,io/Pn,oo = Pio,io/pii,oo remains un- 
changed. Note that the above analysis was based on cer- 
tain symmetry assumptions so it does not exhaustively 
cover all possible local Gaussian operations. Still, it pro- 
vides a strong indication that local Gaussian operations 
cannot help to reduce value of e. 



IV. TWO-COPY DE-GAUSSIFICATION 

We now present a two-copy de-Gaussification scheme 
that reduces the factor e, see Fig. 2. The initial state 
shared by Alice and Bob can be written as Pin,AB®Pin.CD 
with modes A and C belonging to Alice and modes B and 
D belonging to Bob. Both Alice and Bob send their two 
modes through a Mach-Zehnder interferometer formed by 
two balanced beam splitters BSi and BS3 (BS2 and BS4). 
A non-Gaussian operation is applied to each mode inside 
the interferometer. We choose a Fock-state filter Zi = 
n — 1 that eliminates the single-photon state, ^i|l) = 0. 
Finally, the output modes C and D are projected onto a 
superposition of vacuum and single-photon state 



1 



:(9|0)+|1)), 



(20) 



(C + cosh(2s))2-S'2' 



where we assume real q. By changing q we can con- 
trol the amount of entanglement of the final distilled 
Gaussian state. Projection onto state ((20|) can be ac- 
complished by a sequence of a coherent displacement 
D{q), single-photon subtraction, another coherent dis- 
placement D^{q) and projection onto vacuum (29| . In 
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mathematical terms, we have 

{q\^{0\D^{q)aDiq)^{0\ia + q). (21) 

To grasp the principle of the scheme shown in Fig. 2 
let us first consider distillation of truncated two-mode 
squeezed states 



|A> 



1 



(|00)+A|11)) 



(22) 



transmitted from the central source to Alice and Bob over 
two identical lossy channels with intensity transmittance 
T. In the basis |00), |10), |01), |11), the density matrix 
of state shared by Alice and Bob reads 



Pin = 




(23) 



where R = I — T. Parameter e of a Gaussian state that 
would be obtained from this resource by a standard Gaus- 
sification procedure as shown in Fig. 1 can be directly 
calculated using Eq. ([TS]) and we obtain e = (1 — T)A, 
as expected. Suppose now that two copies of the state 
([23| are used as an input of the de-Gaussification scheme 
of Fig. 2. It turns out that the interference on balanced 
beam splitters followed by filters h—1 suppresses all con- 
tributions stemming from the loss of an odd number of 
photons (one or three) during the transmission. After 
some algebra we obtain density matrix of the state after 
de-Gaussification, 



Pab 



^4(1 + 2A2i?2) + A4(g2i?2 _^ 2^2)2 



g4(l + A2i?2)2 









X^T^R^q^ 




\^l^q'- 





(24) 



We can determine the parameter e' of this state similarly 
as before for p-m and we find 

e' = (l-r)2A2. (25) 

This means that e' — and since e < 1 we have e' < e 
as desired. 

We next prove that this property remains true even 
for input symmetric Gaussian states ([T]). The overall 
transformation induced by the filtration can be expressed 

as 

P'aB = {F\^' Fj^D)Pi^.AB Pin^CoiFAC ® Fbd), (26) 

where the local two-mode filters are identical, Fac — 
Fed, and e.g. the filter on Alice's side can be expressed 
as 

Fac = uls.AcifiA - 1)(% - l)%s,Ac|g)c. (27) 

The balanced beam splitter transforms input annihilation 
operators into their balanced combinations, 

^BS,ACa^BS,AC = -^(a-fc), 

uI,s.ac^Ubs,ac = :^('^~^)' (28) 

With the help of transformation rules (pS)) we obtain after 
some algebra 

Fac = |(nA-l)(nA-4)|0)c + |nA(nA-5)|l)c 



--a2(gV2|2)c + V6|3)c) . (29) 

The filtration effectively replaces projection onto a Fock 
state \n)A with projection onto state F\n)A- In particu- 
lar, for the three lowest Fock states we obtain the map- 
ping 

\Q)a ^ 9|0,0)ac, 
^ -|1,1)ac, 

|2)a ^ -i|2)^(q|0)c + 3|l)c) 

-i|0)^(q|2)c + V3|3)c). 

Similar formulas hold for filter on Bob's side. Due to 
the symmetric structure of the Gaussian state (dJ, the 
only nonzero density matrix elements of state p' relevant 
for determination of the asymptotic Gaussian state after 
Gaussification procedure are the following, 

/ 4 2 

Poo, 00 ~ 9 Poo, 00 1 

P'io,io ~ Poi,oi ~ 9^Pio,iO' 

Pn,oo = Poo,ii = Q^Pii.oo- (30) 

If we insert matrix elements (|30p into Eq. we find 

that indeed the parameter eout of a Gaussian state ob- 
tained by Gaussification of p' is equal to the square of ein 
of the initial Gaussian state pin, 

Cout = efn- (31) 
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FIG. 3: Parameter e'^^ of distilled Gaussian state after A*' 
stages of nested entanglement distillation protocol is plotted 
as a function of ein. The numerical labels indicate number of 
stages N of the protocol. 



By iterating the whole entanglement distillation pro- 
tocol several times the parameter e can be reduced to an 
arbitrarily small value. Such a nested distillation pro- 
tocol is similar to a quantum repeater scheme (sol . l3ll | . 
Alice and Bob divide the shared copies of states pin into 
many blocks. They perform distillation on each block, 
thereby obtaining several copies of distilled state /5°°. 
These states are then used as an input for the next stage 
of the nested distillation scheme. After N stages of this 
protocol the parameter e is reduced to 



V. CONCLUSIONS 

In this paper we have proposed an entanglement dis- 
tillation scheme for symmetric Gaussian states that al- 
lows to asymptotically extract a pure entangled Gaus- 
sian state from an arbitrary input entangled symmetric 
Gaussian state. The present scheme is a modified and ex- 
tended version of the entanglement distillation protocol 
for Gaussian states developed by Browne and cowork- 
ers m, [13]. We employ a two-copy de-Gaussification 
procedure and a nested entanglement distillation scheme 
where outputs of a full iterative distillation protocol are 
used as inputs of the next stage of the nested scheme. To 
make the presentation as transparent and comprehensi- 
ble as possible we have considered an idealized scenario 
where all local operations and measurements are imple- 
mented perfectly, i.e. without any errors. In practice, 
both photon addition and subtraction as well as inter- 
ference on beam splitters and projection onto vacuum 
state will suffer from some errors. Also, any practical 
implementation of Gaussification will be only approxi- 
mate, with finite number of iterations. The number of 
iterations will influence both the convergence and suc- 
cess rate of the protocol. An analysis of the influence of 
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This dependence is plotted in Fig. 3. We can see that 
the protocol quickly converges to e = corresponding to 
effective channel transmittance Toff — 1 and a pure state. 

Symmetric two-mode Gaussian states can be charac- 
terized by their purity V and entanglement of formation 
Ef that can be calculated analytically for this class of 
states [13]. In Fig. 4 we plot purity and entanglement of 
formation of Gaussian state obtained after stages of 
the nested distillation protocol as a function of the initial 
effective channel transmittance T. At each stage of the 
protocol, q is adjusted such that the effective two-mode 
squeezing given by Eq. remains constant, r = 1. Each 
stage of the protocol thus increases the entanglement of 
the state, see Fig. 4(b). Purity of the state is also in- 
creased, except for the region of very low T, where state 
purity can be actually reduced after the first stage of the 
protocol. Nevertheless, after subsequent stages purity 
increases and approaches asymptotic value V = \. 




0.0 0.2 0.4 0.6 0.8 



FIG. 4: Purity (a) and entanglement of formation (b) of dis- 
tilled Gaussian state after stages of nested entanglement 
distillation protocol is plotted as a function of the effective 
channel transmittance T for a fixed effective two-mode squeez- 
ing constant r = \. The numerical labels indicate number of 
stages A'' of the protocol, dashed lines represent values of V 
and E; for the initial state before distillation. 
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all these effects on the performance of the protocol will 
be a subject of the future work. 

The entanglement distillation protocol discused in this 
paper combines several elementary building blocks, each 
of which has been already successfully demonstrated ex- 
perimentally. The Gaussification of two and three copies 
of two-mode entangled states of light by their local inter- 
ference on beam splitters and Gaussian measurements 
on some output ports has been employed to experi- 
mentally distill entanglement of phase-diffused two-mode 
squeezed states [tI l33j| The required non-Gaussian oper- 
ation Z — "fi — 1 — 2a^d — da' is a linear combination of 
two sequences of photon addition and photon subtraction 
that has been recently implemented experimentally for 
the purpose of direct experimental testing of fundamen- 
tal quantum commutation relations (28] . In view of these 
recent achievements, a proof-of-principle demonstration 



of the entanglement distillation scheme proposed in this 
paper appears to be within the current experimental lim- 
its. However, the purely optical scheme is likely to exhibit 
exponential decrease of success probability with grow- 
ing number of distilled copies. This unfavourable scal- 
ing could be avoided with the help of a quantum mem- 
ory [s^ . Issj , similarly to the quantum repeater schemes 
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